Dynamical noise and avalanches in quasi-static plastic flow of amorphous solids 
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We build a mean-field model of plasticity of amorphous solids, based on the dynamics of an 
ensemble shear transformation zones, interacting via intrinsic dynamical noise generated by the 
zone flips themselves. We compare the quasi-static, steady-state properties for two types of noise 
spectrum: (G) Gaussian; (E) broad distribution derived from quadrupolar elastic interactions. We 
find that the plastic flow proceeds via avalanches whose scaling properties with system size are 
highly sensitive to noise tails. Comparison with available data suggests that non-affine strain fields 
might be of paramount importance in the small systems accessible to molecular simulations. 



Decades of efforts to explain plasticity of amorphous 
solids, and more generally of jammed disordered sys- 
tems (foams, granular media, colloidal glasses.. .), have 
converged towards general agreement on the nature of 
the elementary dissipative events in these highly mul- 
tistable systems. They involve sudden rearrangements 
of small clusters comprising a few basic structural units 
(molecules, bubbles, grains,...). Each 'flip' can be 
viewed as the shear-induced transformation of an ellip- 
soidal Eshelby inclusion and generates a long-range elas- 
tic displacement field with quadrupolar symmetry. In 
spite of this progress, relating the macroscopic stress re- 
sponse to the statistics of elementary events remains a 
challenging question. 

Recently, twophenomenologies of plasticity of jammed 
systems (SGR H and STZ 2]) have been proposed. For 
zero-temperature systems, they derive constitutive equa- 
tions based on a few common assumptions for the shear- 
induced dynamics of an ensemble of non-interacting 
shear 'zones' (or 'traps'). Namely, the flippable zones 
are advected by the external drive at shear rate e to- 
ward their instability threshold; meanwhile, their elastic 
energy fluctuates under the effect of a mechanical noise, 
assumed to be described by a fixed effective temperature. 
Flips can thus be triggered before zones reach the insta- 
bility threshold, with strain-rate independent^ Arrhenius- 
like jump probabilities. In this sense noise, in these mod- 
els, can be viewed as extrinsic. 

Given that noise has a crucial bearing on the form 
of the resulting constitutive laws, it would be desirable, 
however, to relate phenomenological models to micro- 
scopic arguments. Following Argon et al 0, consider 
a zone far from its threshold. As time elapses, its strain 
is advected by the external drive. In steady state, flips 
occurring at random sites elsewhere in the material emit- 
at an average rate ~ e-elastic signals which propagate 
at sound speed and upon arrival, shift the local strain. 
These fluctuating shifts, originating from spatially ran- 
dom sources, constitute the dynamical ^ intrinsic^ me- 
chanical noise. In view of the long range of elastic cou- 
plings, one can expect it to give rise, as in comparable 



systems with pinning and long range interactions to 
avalanches, i.e. 'instantaneous' series of events triggered 
by a single initial flip. The existence of correlated flips 
may have a bearing on the macroscopic plastic proper- 
ties, including the emergence of shear localization. 

In this letter, we analyze a mean-field model for the 
quasi-static dynamics of a random 2D ensemble of N 
zones, driven by shear and interacting only via intrinsic, 
dynamical noise. We first clarify the limits of the quasi- 
static regime in which flips can be assumed instantaneous 
and avalanches are well-separated. We then construct 
two versions of the model, differing in their assumptions 
about the noise spectrum: model G with Gaussian noise 
and model E more realistically based upon the quadrupo- 
lar elastic field In both cases, we find, in agree- 
ment withprevious numerical simulations of molecular 
glasses 0, B 01 that flips cluster into avalanches of 
average size (n) increasing as a power a of system size 
N. As expected, the average stress drop decreases with 
N ensuring proper convergence towards a smooth stress 
curve. The two noise spectra, however, lead to drastic 
differences in several important features, including: the 
average shear stress, the distribution of avalanche sizes 
and its scaling properties, and the scaling exponent a. 
For model G, simple arguments based upon an approx- 
imate master equation provide a reasonable prediction 
for a. A similar argument fails for model E, probably 
due to the presence of tails larger that Gaussian, even 
though all the moments of the elastic noise are finite. 
Confrontation with molecular simulation results points 
towards the need for further developments, including in 
particular non-affine strain field effects. 

We consider an ensemble of A'' identical zones of size a, 
randomly distributed with a fixed density a^/d^, in a 2D 
elastic medium of lateral size L. The elastic state of each 
zone is characterized by an internal strain , which mea- 
sures the departure from its zero stress state. The e^'s 
lie below a common stability threshold ec. The system is 
driven by external shear at rate e, which advects all e^. 
When a zone reaches Ec, it disappears and another one 
is created at an uncorrelated position, with zero initial 
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stress (hence zero local strain). The duration of a flip, 
controlled by acoustic radiative damping, is to a/cg 
with Cs a sound speed. A flip at site emits an acoustic 
signal which propagates throughout, and modifies the lo- 
cal strains of all other zones j, which adjust over a time 
~ ^0, to this space-dependent shift 6e{rij). Following Pi- 
card et al 13 , this elastic field has quadrupolar symmetry, 
hence zero average: it constitutes a noise acting on the 

In steady state, the average flip rate in the whole 
system is i?flip — 5t^l^ — Ne/tc- In the absence of 
avalanches, the noise correlation time is and the quasi- 
static regime corresponds to 5ifiip to i-e. e <C eflip — 
-^j^- For a glass-like system, with Ec ~ 1%, a ~ Inm, 
zone density a? /(P ^ 10~^, lateral size L ~ Imm, this 
yields the loose criterion e ^ ls~^. Now, an elastic noise 
signal may drive some e^ 's beyond Ec, hence trigger sec- 
ondary flips, thus initiating an avalanche, whose duration 
ta is set by sound propagation. For a conservative esti- 
mate we take the average distance between successive 
flips to be L. This leads to ta = {n)L/cs with (n), the 
average avalanche size. The average avalanche frequency 
is thus i?a = ^ = (n)l ' quasi-static condition 

becomes 11]: e ^ = ^ eflip. Avalanches drastically 
reduce the quasi-static range. A second important point, 
usually ignored, is the rapid shrinking of this range with 
system size. 

Models: In the quasi-static regime the dynamics re- 
duces to steady drift of all ei's at fixed e, interspersed 
with instantaneous avalanches. We can eliminate time 
and define our models following a two-step algorithm: 
(i) We start from a configuration where all e,; < Ec- We 
identify em = Max(ei) and shift the macroscopic strain, 
hence all e^'s, by Aedrift — — ^m- (ii) Zone m flips, 
i.e. is removed while another one is introduced at zero 
strain. The resulting elastic signal is modelled by random 
shifts Sci of all other e^'s. Here, we make the mean- field 
assumption that the ^e^'s are independent random vari- 
ables, which amounts to ignoring the correlations due to 
the spatial zone arrangement. The elastic nature of dy- 
namical noise then enters via the spectrum of ^e^'s (see 
below). This first flip yields a new configuration e|^'. If 
em^ < Cc we are back to step (i). Otherwise, an avalanche 
starts: we count the number z ^ qi oi zones that flip 
at this stage and treat their noise signals successively. 
If the first signal triggers q2 new flips, z is updated as 
z — !■ z + 52 ^ 1, etc. The avalanche stops when z — 0; 
its size n — 1 + '^qi. We implement this algorithm with 
two specifications of dynamical noise. 

Model E is based upon the results of Picard et 
al 13 : the elastic strain at position (r, 0) from a flip 
is: f ° cos(46'), where 9 is measured from a prin- 
cipal axis of the macroscopic strain, and the elemen- 
tary strain Aep/ec = tt L'^ /2N a'^ . Each instance of the 
noise dei is obtained by uniform Monte-Carlo sampling of 



(r, 9) L[x] — TT, 7r[, where d is the average distance be- 
tween nearest zones and = N d'^. Note the size depen- 
dence of this elastic noise, of variance A^v = '^A.ge'^/N, 
where, in our circular geometry, the factor Aq = I/tt^. 
One easily checks that the noise distribution p^iSe) ex- 
hibits a l/(5e^ behavior up to the cut-off 2a^Aeo/7r(i^ 
beyond which it vanishes, so that all its moments are 
finite. 

Model G assumes that noise is simply Gaussian with 
the same variance A^r. It will provide a test of the influ- 
ence of tail shapes in p. 

Numerical results: In order to study the statistical 
properties of such models in steady state, we elimi- 
nate the initial transient (e < 2 Ec) and perform er- 
godic averaging over long time intervals involving > 10^ 
avalanches. We plot on Figure ^ the variations of the 
average avalanche size (n) with system size N . The data 
for model E fit closely the power law (n)^ ^ N"'^, with 
aE = 0.147, over the whole range (two decades). The 
data for model G are consistent with an asymptotic (see 
insert) power law behavior (n)^ ~ N"'^ , with the trivial 
exponent ac = 1/2. 

We also compute the ergodic average of the macro- 
scopic stress cr. Since the macroscopic strain equals the 
average zone strain, a — 2 p (ei) (with p the shear modu- 
lus). We find for ^ = a/ {2 pec) the values /3g = 0.40 
and f3E = 0.42. The 5% difference, though small, is 
numerically significant. We also check, in agreement 
with the above power law scalings, that stress fluctu- 
ations decrease rapidly with system size, the ranges of 
the noisy stress curves becoming completely separated 
for N > 32000. Different (3 values imply differences in 
the stationary distribution p of e^. p{e) is represented 
on Figure |2] for both models and various values of N. 
For each model, p converges rapidly almost everywhere 
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FIG. 1: Average avalanche size (n) vs TV for models G (trian- 
gles) and E (circles). The dashed line has slope 1/2. Insert: 
{n}N-"/{n){NQ)No°' vs N/No with A^o = 250. 
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the evolution of p(e) by the approximate master equation: 



FIG. 2; Steady state distribution p of zone strains vs e/tc 
for: (a) model E; (b) model G. Inserts: blow-ups of the near- 
threshold region. The arrows indicate increasing values of 
iV = 250 X 2" = 0, .. .,7). 



towards a limit curve, which explains that the /3 values 
vary by less than 0.4% (model G) and 0.02% (model E) 
when N increases from 1000 to 32000. While pc and 
Pe are similar in most of the e range, they present sig- 
nificant differences in two regions, e ^ and e ~ Cc- 
The peak 'm pe results from refeeding zones at e = 
after flips. The larger an avalanche, the more the corre- 
sponding peak is broadened by ulterior flips within the 
avalanche itself. We therefore attribute the washing out 
of the peak for model G to the fact that it exhibits much 
larger avalanches. More significant for our purpose is 
the detailed behavior of p near threshold, which reflects 
avalanche statistics. In particular, pc = p{^c) is directly 
related to (n). Indeed, the average flip and avalanche 
rates verify i?fiip — i?a {n)- On the other hand, the ad- 
vected zone flux in steady state reads /adv. — £PcN. It 
comprises all zones initiating avalanches (while the "de- 
scendants" of each of these mother zones feed the noise 
contribution to the total flux i?flip). Hence /adv. — Ra, 
from which pc — {{n) ec)~^. Since the data indicate that 
(n) diverges for N ^ oo, they also indicate that the ab- 
sorbing boundary condition pc = should hold asymp- 
totically here. This relation provides a consistency test 
of our calculations. We determine Pc with the help of 
a second order polynomial extrapolation near Ec, with 
sampling intervals Ae/ec = 1.25 x 10^^. We flnd that 
the relation Pc{n) Cc = 1 holds within 1% for model G, 
and 3% for model E. 

Discussion In order to try and obtain analytical esti- 
mates for the above avalanche size scalings, we describe 
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where w is the single flip transition probability: w{6e) = 
N e p(5t) / Cc with p the above-defined noise distribution. 
Here F = ^_^d5tw{5e). The delta term, propor- 
tional to the normalized zone flux f /N ~ e/ec, accounts 
for post-flip reinjection and ensures the conservation of 
zone number. In this approximation, advection oper- 
ates between all single flips, which amounts to neglect- 
ing intra-avalanche time correlations. This leads to an 
average advective e-shift during an avalanche Acadv. ^ 
e/-Ra ^ Ec {n)/N, to be compared with the average dif- 
fusive broadening Aediff: ^ ^/(ri) Ajv with Ajv '^-^ e^/N 
the noise variance. Hence, Aeadv./Aedig. ~ \/ (n) /N ^ 
jy-{i-a)/2 ^]^[(.h suggests that our approximation should 
improve in the large N limit. Integration of equation 
in steady state yields: 



dep(e) / de'w(e — e') 
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T7 = — = 
N er 



Since w is peaked around zero, we expand p{e) close to 
Ec to first order: p(e) ^ Pc + — ^cjV'c- Using pc — 
((n) ec)~^, we obtain for the average avalanche size: 



(n) 



Np'^ ((5e 



n -1 



(2) 



where the (semi)-moments (<5e'')-|_ = d((5e) w((5e) (5e'". 
For both model, ((5e^ 



and 



/(oo) 



= Aoel/N, while ((5e)f^ = 
If p'^ converges to- 



Scc log N 



N 



equation ^ predicts that, for 



wards a finite value pc 
large systems, 

- for model G: (n) - N^^^ 

- for model E: (n) ~ log TV 

While this prediction accounts satisfactorily for the nu- 
merical results for model G, we have checked (see also 
Fig. nj that the log scahng is ruled out by our data. The 
reason for this failure is illustrated by the insert of Fig.[21- 
(a). For model E, we find that, for increasing N, p be- 
comes increasingly steep in the near vicinity of ec- In the 
TV-range investigated, we see a marked, non-saturating, 
increase of \p'J, suggesting a possible divergence, higher 
derivatives increasing even faster. This highly singular 
behavior clearly invalidates the above truncated expan- 
sion of p for model E. Conversely, for model G, we find 
numerically that p'^ Aq does exhibits convergence to- 
wards ~ —1.3. 

This discussion points towards our central result: the 
slowly space decrease typical of elastic interactions cru- 
cially affects avalanche behavior. This holds not only 
for their average size, but, as well for their size distribu- 
tion n7(n). As appears on Fig. |31 for both models, w(n) 
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FIG. 3: Distribution uj{n) of avalanche sizes, for model G 
(left) and E(riglit), and for all system sizes. Inserts: (n) vj{n) 
vs n/{n). The arrows indicate the direction of increasing sizes. 



sue of localization, full simulations of model E will be 
necessary. 

The second series of questions arises when confronting 
our results with those of |^ . Maloney et al find in the 
quasi-static regime that (TV) scales approximately as \/N. 
Moreover their avalanche distribution obeys quite well 
the above-mentioned scaling. Analogous results emerge 
from the 3D simulation of an amorphous metallic alloy by 
Bailey et al . This seems to plead in favor of Gaussian 
noise. This a priori surprising conclusion raises a further 
question. Recently, Leonforte et have shown that 

the elastic response of amorphous solids self-averages into 
the continuum elastic response only beyond a length scale 
^ of order ~ 20 molecular diameters Oq. For r < ^, the 
elastic response is dominated by non-affine effects. So, for 
the existing molecular simulations focusing on avalanche 
dynamics, where L is limited to < 50 ao units, noise tails 
are very likely to be controlled by elastic non-affinity. 
Incorporating such effects into models of intrinsic noise 
thus appears as an important, though intricate, issue. 



decreases exponentially for n 3> (n). However, a fun- 
damental difference appears when attempting to rescale 
w{n) as ro(n) = (l/(n)) /(n/(n)) (inserts of Fig. |3J). For 
model G, the quality of the collapse, good for n > (n) 
at all N, improves with system size in the small n/{n) 
range. On the contrary no such scaling holds for model 
E. 

These results on dynamical noise effects, which remain 
preliminary, indicate two routes for further investiga- 
tions. 

Our mean-field approximation wipes out from the start 
the directional effects arising from the quadrupolar struc- 
ture of the elementary events, which are responsible for 
the preferential avalanche orientations observed in the 
LJ glass simulations by Maloney et al ^ ^ and Tan- 
guy et al 9]. In order to evaluate the robustness of the 
mean-field scalings and also to start addressing the is- 
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